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Abstract: Many early-type stars are in binary systems. A number of them shows radio emissivity with pe- 
riodic variability This variability is associated with non-thermal synchrotron radiation emitted by relativistic 
electrons. The strong shocks necessary to accelerate the electrons up to high energies are produced by the 
collision of the radiatively-driven stellar winds. A study of the non-thermal emission is necessary in order to 
investigate O-star colliding wind binaries. Here preliminary results of our modeling of the colliding winds in 
Cyg 0B2 No. 9 are presented. 

1 Introduction 

Many OB stars show non-thermal radio emission together with the thermal free-free radiation in the 
same band. The non-thermal emissivity is thought to be due to synchrotron radiation emitted by 
relativistic electrons spiraling around the dipole magnetic field lines of the star. The electrons are 
accelerated up to high energies by strong shocks. 

In recent years a number of early-type stars radiating non-thermal radio emission have been dis- 
covered to be binary systems. This discovery supports the hypothesis of colliding winds from the 
primary and the secondary as the origin of the strong shocks (see Eichler & Usov 1993, Dougherty 
et al. 2003, Pittard et al. 2006). The observed synchrotron emissivity is characterized by a periodic 
variability directly connected with the orbital phase of the binary system. 

One of these binary systems is Cyg 0B2 No. 9, which was detected as a non-thermal emitter in 
1984 (see Abbott et al. 1984). By studying the observed VLA (Very Large Array) radio fluxes at 
3.6, 6 and 20 cm (Van Loo et al. 2008) from Cyg 0B2 No. 9, Van Loo and collaborators found a 
variability with a period of 2.35 yr. This periodicity confirmed the presence of radio non-thermal 
emission suggesting also a binarity of this system. At the same time Naze et al. (2008) discovered its 
binarity after a long-term spectroscopic monitoring. The observed 6 cm radio flux is shown in Fig.[Tl 
The orbital parameters of Cyg 0B2 No. 9 were derived for the first time by Naze et al. (2010). 

The study of the synchrotron emission at different orbital positions is relevant to constrain stellar, 
orbital and wind parameters, such as the mass loss rates from the stars and the clumping/porosity, 
and in general to investigate the nature of the stellar winds. As a matter of fact the mass loss rate 
determination in massive stars is presently one of the main problems. This problem has become 
relevant due to the clumping/porosity of the stellar winds which affects the estimation of the mass 
loss rate. To determine the emitted synchrotron radiation and its absorption by the stellar wind plasma 
permits to estimate the amount of porosity in the wind. 

■ "The multi- wavelength view of hot, massive stars"; 39^^ Liege Int. Astroph. Coll., 12-16 July 2010 ■ 



In order to better understand the physics of the non-thermal emission and to constrain the model 
parameters we developed a numerical code. Here the simulated model for Cyg 0B2 No. 9 is pre- 
sented. The theoretical results are compared with the corresponding observations. 
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Figure 1: Cyg 0B2 No. 9. Observed VLA 6 cm radio emission (from Van Loo et al. 2008). 



2 Modeling 

In this section the theoretical model used to obtain the simulated results is described. A first version 
of this model was developed by R.Blomme. The details of the model can be found in Blomme et al. 
(2010). 

The two winds are considered to collide at the contact discontinuity. Its position is assumed to be 
coincident with the two shocks and it is calculated using the Anthokin et al. (2004) equations. The 
electrons accelerated at the shock are followed as they advect away and cool down along the post- 
shock streamlines. We stop following the electrons when their momentum falls below a minimum 
value (see Van Loo et al. 2004), or when they leave the simulation volume. Adiabatic and inverse 
Compton losses are both taken in account. The momenta follow a modified power-law distribution. 

The synchrotron emissivity from the relativistic electrons is calculated along the post-shock stream- 
lines in the orbital plane. The Razin effect is also included. The formula of the emissivity ju{r) is 
reported below as in Van Loo et al. (2005): 

Jv(r) = — / dpN{p,r) —2nsme^^Bf{u,p)smeF[— — — - (1) 
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where u is the frequency, r is the radial position, pq is the minimum momentum, pc is the maximum 
one, p is the particle momentum, N{p, r) is the modified power-law distribution of the electrons, 9 is 
the pitch angle, e and me are the electron charge and mass, c is the velocity of light, B is the magnetic 
field at a distance r, z/g is the critical frequency, and /(z^,p) takes into account the Razin effect. 

The synchrotron emission is calculated in the orbital plane. The third dimension is recovered by 
rotating the orbital plane along the line which connects the two stars. The free-free opacity and emis- 
sion due to the ionized wind material are also included using the Wright & Barlow (1975) equations. 
The theoretical fluxes are obtained at different orbital phases using Adam's method (Adam 1990) that 
solves the radiative transfer equation in a very simple way. 

The orbital parameters are provided by Naze et al. (2010) observations, the stellar parameters by 
Martins et al. (2005) theoretical models and the wind parameters by Vink et al. (2001) theoretical 



Table 1: Stellar, orbital and wind parameters for Cyg 0B2 No. 9. 





primary star secondary star 


P (days) 


860.2 ± 5.5 


e 


0.744 ± 0.030 


i (deg) 


48.59 ± 7.45 


u (deg) 


-164.4 ±4.1 


d(pc) 


1820 lb 200 


vq (km/s) 


0.1 




50.72 44.10 


« (-Rsun) 


797.34 ± 203.95 932.40 ± 237. 17 


R {Rsun) 


19.45 19.95 


L (erg s~^) 


2.85-1039 2.48-1039 


Teff (K) 


38612 36801 


f oo (km/s) 


2050 1850 


M (Msun/yr) 


5.79d-06 4.86d-06 



Note. P and d are the orbital period and the distance from the Earth (Van Loo et al. 2008), e and uj are the 
eccentricity and the angle of the periastron (Naze et a. 2010), M^,, R, L, Tcff are the mass, the radius, the 
luminosity, and effective temperature for the star (Martins et al. 2005), i is the inclination angle (obtained 
from a comparison between the theoretical mass A/^ calculated in Martins et al. 2005 and the projected 
observed one Mi, aini^ from Naze 2010), a is the semi-major axis of the orbit, vq is the velocity of the 
wind at the surface of the star (Blomme et al. 2010), Voo and M are the velocity of the wind at large 
distance from the star and the mass loss rate (Vink et al. 2001). 



models. We choose spectral type 05I-I-06I: the Martins' stellar parameters for this spectral type are 
consistent with the spectroscopic observations and the chosen spectral type is in agreement with the 
05-1-06-7 spectral types suggested by Naze et al. (2008). Current spectroscopic information does not 
allow a reliable determination of the luminosity class. 

The knowledge of the spectral type is necessary to select the values of parameters that cannot be 
obtained by the observations. These parameters are used as input to our model. The chosen values of 
the parameters are reported in Tabled] 

3 Theoretical results: the emissivity and model light curve at 
6 cm 

The 6 cm emissivity images at periastron and apastron are shown in Fig.|2l The 6 cm emissivity at pe- 
riastron is much more concentrated around the contact discontinuity than at apastron mainly because 
the Razin effect is more important in the periastron high-density region. The synchrotron radiation is 
however not negligible in both cases due to the extension of the emitted area. This extended region 
explains the synchrotron emission that was observed in OB stars systems during the last years. 

In Fig. |3] the preliminary 6 cm light curve from our modeling is presented. The variability in the 
radio flux is clearly visible as observed with VLA. This variability linked to the orbital period is the 
fingerprint of the non-thermal radiation. Compared to the observations (see Fig.[T]), the theoretical 
fluxes are much too high and the maximum occurs too early. This overestimate of the flux could 
be produced by too high a number of particles being taken into account. In order to reduce the 
theoretical fluxes we could decrease the fraction of energy transferred from the shock to the relativistic 
particles. In the present work the assumed value was the typical 0.05 accepted for a strong shock 
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Figure 2: Cyg 0B2 No. 9 simulated results. Left panel: the 6 cm emissivity at periastron. Middle 
panel: the eccentric orbit of the binary system. Right panel: the 6 cm emissivity at apastron. Along 
the axis the distances from the apex are in unit of Rsun- The color bar represents the emissivity 
in erg/cm^ s Hz with a logarithmic scale. The arrows link respectively the periastron and apastron 
positions with the corresponding emissivity pictures. 

following Eichler & Usov (1993). Another option could be to increase the index n of the p^" particle 
momentum distribution function and thus to decrease the shock strength, or even to diminish the star 
surface magnetic fields that we assumed to be 100 Gauss. Before tuning these model parameters it is 
however more urgent to evolve directly the hydrodynamical equations in order to correctly position 
the shocks and obtain a realistic value for their strengths, and to take into account the instabilities and 
inhomogeneities that influence both the thermal and non-thermal emission. Then it would be relevant 
to include the orbital motion. 

Further investigation and changes in the model will be thus necessary to improve quantitatively 
the simulated results. 
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Figure 3: Cyg 0B2 No. 9. Simulated 6 cm light curve. 



4 Future work 



Future work will consist first of all in solving directly the hydrodynamical equations and including 
the orbital motion. Then it will be necessary to explore the parameter space of our model to obtain a 
better agreement with the corresponding observations (included the stellar spectral types). 
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